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Abstract 

We present the results of a numerical calculation of semi-leptonic form 
factors relevant for heavy flavour meson decays into light mesons. The 
results have been obtained by studying two- and three-point correlation 
functions at (3 = 6.4 on a 24 3 x 60 lattice, using the Wilson action in the 
quenched approximation. 

From the study of the matrix element < K~\J il \D Q > we obtain f^(0) = 
0.65 ± 0.18. From the matrix element < K*°\J^\D + > we obtain U(0) = 
0.95 ± 0.34, Ai(0) = 0.63 ± 0.14 and A 2 (0) = 0.45 ± 0.33. We also obtain 
Ai{q 2 max ) = 0.62 ± 0.09, V(0)/Al(0) = 1.5 ± 0.28 and A 2 (0)/Al(0) = 
0.7 ± 0.4. The results for /£(0), V(0) and A x (0) are consistent with the 
experimental data and with previous lattice determinations with larger 
lattice spacings. In the case of ^(0) the errors are too large to draw any 
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firm conclusion. We also show, with the help of the heavy quark effective 
theory (HQET), that it is possible to extrapolate the form factors to the 
B meson. Our calculations show that the form factors follow a behaviour 
compatible with the predictions from HQET. Within large uncertainties, 
our results suggest that AijA\ increases with the heavy quark mass. For B 
mesons A2/A1 can be as large as 1.5-2.0. We also get very rough estimates 
for the partial decay widths T(B -> ttIvi) = \V ub \ 2 (l2 ± 8) x 10 12 s _1 and 
T(B — ► plv{) = |Kb| 2 (13 ± 12) x 10 12 s _1 , which can be used to give upper 
bounds on the rates. 



1 Introduction 



Semi-leptonic decays of heavy-light mesons have attracted considerable interest 
in the past years as they play a crucial role in the determination of the Cabibbo- 
Kobayashi-Maskawa mixing matrix and in the understanding of weak decays. 
Moreover, the study of the dependence of the form factors on the heavy quark 
mass checks the validity of the scaling laws, predicted by the heavy quark effective 
theory, in the range of masses corresponding to D and B mesons. 

There is increasing evidence, that lattice QCD allows ab initio quantitative 
predictions of weak decay matrix element sjlj. Exclusive semi-leptonic decay chan- 
nels of heavy flavour mesons have been studied in a series of papers at a 
lattice resolution of a -1 ~ 2 — 3 GeV. In order to improve the control over dis- 
cretization errors, we present here a study on pseudoscalar-pseudoscalar as well 
as pseudoscalar-vector weak form factors, on a lattice 24 3 x 60, at /3 = 6.4, which 
corresponds to a -1 ~ 3.6 — 3.7 GeV. This is a continuation of our recent work|| 
on the leptonic decay constants. 

The main results of the present study are given in the abstract and in tables |§ 



T0| and O. These results have been obtained by extrapolating the matrix elements 



in momentum transfer, in the light and heavy quark masses. We have included in 
tables [TJ-|S] all the numbers obtained from the direct study of the two- and three- 
point correlation functions, computed at several values of K w . These tables may 
be useful to check our results. We also give the details of the extrapolation, in 
order to allow the reader to reproduce the form factors in table from the meson 



masses and the form factors of tables In section |3.2| we show that it is possible 
to extrapolate the form factors to the B meson using the Heavy Quark Effective 
Theory (HQET). In this way one can get also informations on the corrections to 
the infinite mass limit. Due to large statistical errors, we are only able to give a 
rough estimate for the branching ratios of the processes B — > 7r and B — > p. 



2 Description of the Calculation 

2.1 Lattice Setup 

We work with the standard Wilson action at /3 — 6.4 for the gauge fields and 
the quark propagators 0, in the quenched approximation. We have generated 15 
independent gauge field configurations on a 24 3 x 30 lattice, separated by 500- 
1600 sweeps, using the overrelaxed algorithm |l(|. The 15 configurations were 
produced in groups of 5, with three independent initial conditions. In each of 
the three cases the first of the useful configurations was obtained after an initial 
thermalization of at least 3000 sweeps ( 500 with the Metropolis algorithm and 
2500 with the overrelaxed algorithm ). 

On each configuration we have computed the quark propagators for 7 different 
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values of the Wilson hopping parameter Kw, corresponding to "heavy" quarks, 
K H = 0.1275, 0.1325, 0.1375, 0.1425, and "light" quarks, K L = 0.1485, 0.1490 
and 0.1495. Periodic boundary conditions on a 24 3 x 60 lattice have been imposed 
in the calculation of the quark propagators, by using appropriate combinations of 
periodic and antiperiodic (in the time direction) quark propagators calculated on 
a 24 3 x 30 lattice [|ll|]. Further details on the lattice calibration, fitting procedures, 



mass spectrum, extraction of matrix elements of local operators between the 
vacuum and meson states, e.g. < M P \Qj 5 q\0 >, can be found in ref.||. 

2.2 Form Factor Evaluation 

From the study of three-point correlation functions 0-0, one extracts the weak 
current matrix elements for a given momentum transfer: 

< K\J,\D >= ( PD+PK - M °- M K q ) f+tf) + Ml- M% gM ^ (1) 

r 2V(a 2 ) 

<K* r \J,\D> = 4 ™^ e^plp s K * + i(M D + m^*) A(g 2 )^ 

- h^ p "+^™*^> (2) 

where q is the momentum transfer, q = po — Pk or q = po — Pk*, P = Pd + Pk* 
and e@ is the polarization vector of the K*. fp°, V, and A are dimensionless 
form factors in the helicity basis, see for example 0. From the matrix elements 
([I]) and (fj), by varying the Lorentz component of the current, the meson momenta 
and the K* polarization, one can extract the form factors. 

In the following we give those details of the calculation which are specific to 
semi-leptonic decays and cannot be found in ref . || . 

The matrix elements have been computed for a D meson at rest using a 
pseudoscalar density as source, at a time distance (to — tx,K*)/a = 28. The 
position of the light meson source is fixed in the origin and we have varied the 
time position (in the interval tj/a = 12 — 16) and the momentum of the weak 
current. For the iT-source we have also used the pseudoscalar density and for 
the K* the local vector current. For the weak axial current we have chosen the 
local operator Z J 4<5(x)7 At 7 5 g(a;), while for the weak vector current we have used 
both the local (Z v Q(x)'j^q(x)) and the "conserved " currents^. Za and Z v are 
the renormalization constants of the axial and vector current respectively fll2|, |i~3|j . 



The values quoted for the form factors are based on the choice Za = 0.88 and 
Z v = 0.840-0 which will be justified below. 

The two- and three-point correlation functions have been computed for differ- 
ent spatial momenta which are allowed by the lattice discretization and volume. 

3 The latter would be conserved on the lattice in the limit of degenerate quark masses [D| . 
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K 1 = K 2 


momentum 


E 5 


E 5 


Z 5 x 10 3 




n n n 


n 97 (1 \ 

U.Zii 111 


n 97(1 \ 




0.1485 


10 

-L, <J, U 


0.39(2) 


0.38(1) 


5.7(1.6) 




1,1,0 


0.41(4) 


0.46(1) 


2.5(1.6) 




0,0,0 


0.23(1) 


0.23(1) 


3.3(0.3) 


0.1490 


1,0,0 


0.37(3) 


0.35(1) 


5.4(2.5) 




1,1,0 


0.37(5) 


0.44(1) 


1.7(1.3) 




0,0,0 


0.19(1) 


0.19(1) 


2.9(0.4) 


0.1495 


1,0,0 


0.36(5) 


0.33(1) 


5.3(5.2) 




1,1,0 


0.32(5) 


0.42(1) 


1.1(0.9) 



Table 1: Pseudoscalar meson masses and energies in lattice units and Z 5 from a 
fit to two-point correlations at different momenta and quark masses. The mesons 
are composed of two degenerate light quarks, and the Wilson parameter is given 
in the first column. 



By defining Pk,k* = 2n/ (La) x (n x , n y , n z ), with L = 24, we have used for Pk,k* 
the following assignments (0, 0, 0), (1, 0, 0), (0, 1, 0), (0, 0, 1), (1, 1, 0), (1, 0, 1), 
(2, 0, 0), and (1, 1, 1). We have combined different correlation functions which 
can be related by the spatial orthogonal group and parity. In all we have five 
independent momenta which in the following we will denote by (0, 0, 0), (1, 0, 0), 
(2,0,0), (1,1,0) and (1,1,1). 

Large momenta in lattice units imply large systematic errors. In order to 
monitor lattice artefacts at large momenta, we have fitted the light pseudoscalar 
meson two-point functions to the asymptotic form (at large time distances): 

G 5 (t,p) = ^e^<P 5 (x,t)P 5 t (0,0)> 

X 

-> Z 5 — ^_cosh(P 5 (t-T/2)), (3) 

where P5 = q(x)^z,q'(x), T is the time lattice length and E§ the meson energy. 
The corresponding formula for the vector meson can be easily derived: 

Gn&p) = 5>** < VKx,t)V/(0,0) >^ 

X 

{-9H + If K^^j^ cosh(P y (t - T/2)), (4) 

where Vi = (j(x)'yiq' (x) . 

An estimate of the lattice effects is obtained from a comparison of E$, obtained 
from eq.(^|), with E 5 = JM^+p 2 , where M5 is the value of P5 found from 
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K 1= K 2 


momentum 


Ey 


Ey 


Zyy x 10 3 




n n n 


U.OO^l J 


n n ( 1 ^ 

U.OO^l ) 




0.1485 


10 

-L, <J, U 


0.42(1) 


0.42(1) 


1.6(2) 




1,1,0 


0.51(1) 


0.49(1) 


2.1(5) 




0,0,0 


0.30(1) 


0.30(1) 


1.2(1) 


0.1490 


1,0,0 


0.40(1) 


0.40(1) 


1.3(2) 




1,1,0 


0.50(2) 


0.48(1) 


2.0(6) 




0,0,0 


0.27(1) 


0.27(1) 


0.9(1) 


0.1495 


1,0,0 


0.38(1) 


0.38(1) 


1.1(2) 




1,1,0 


0.51(2) 


0.46(1) 


2.3(8) 



Table 2: Vector meson masses and energies in lattice units and Zyy from a fit 
to two-point correlations at different momenta and quark masses. The mesons 
are composed of two degenerate light quarks, and the Wilson parameter is given 
in the first column. 



a fit at p — 0. Moreover the values of Z 5 obtained by fitting the two-point 
functions at different momenta should be consistent (similarly for vector mesons). 
In analyzing the two-point correlation functions, we found that for (1, 1, 1) and 
(2, 0, 0) the results are too noisy with the present statistics and therefore they 
have not been used. In tables p] and we report, for different light quark masses, 
E 5 (E v ) and Z 5 [Zyy) as obtained using the fit in eq.([3|) (eq.(f|)), as well as E 5 
(Ey). The results reported in tables [l] and § have been obtained by fitting the 
two-point functions in the interval 12 < t/a < 20 Tables [I] and |^ show that for 
the momenta (1,0,0) and (1,1,0) the continuum energy-momentum relation is 
satisfied within large statistical fluctuations. The statistical error increases as 
one approaches to the chiral limit and the agreement appears to be better for the 
vector meson than for the pseudoscalar one. We believe that this is due to our 
limited statistics. 

In order to extract the current matrix elements from the three-point functions 
several approaches may be followed. We have used two different methods denoted 
in the following by "analytic" and "ratio" . To be specific, we only explain these 
methods in the case of < K\ J^D >. The three-point correlation function is given 
by: 

C^ty) = E<°I T 

x,y 

- < + > e-^M-WW (5) 

at large time distances. In eq.(|5|), P§ (Pl D ) is the pseudoscalar density used to 



^(0,0) J, 



^x, t x j 



Pl D {y 



i tyj 



10 > 



Jq-x+ip D -y 
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annihilate (create) a K meson (D meson) and E 5 ' the corresponding energies^. 
The matrix element can be isolated from the ratio: 

R = TO^i x ^ ^ KWD > (6) 

where G^' D are the analog of G 5 defined in eq.(||). 

The two methods to extract the matrix elements are as follows: 

• "analytic": in the denominator of eq.(|B|), we use the analytic expression 
of G§' D given in eq.(|3|), with Z 5 and E 5 (computed from the meson mass) 
taken from the fit to the two-point function at zero momentum. 

• "ratio": we divide the three-point correlation function by the two-point 
correlation functions with appropriate momenta averaged over the same 



configurations. We then multiply the result by the factor \JZ^Z® where 
Z§ and Z® are obtained from the fit at zero momentum to the two-point 
functions corresponding to the D and K propagators. 

With large time distances and a high statistics the two methods should agree, up 
to 0(a) effects. Our limited statistics gives fluctuations in the energy- momentum 
relations, in addition to these 0(a) effects. Therefore the two methods may 
yield slighlty different results and we will take into account the differences in the 
evaluation of the final error, see below. We report in tables [3] and |] the form 
factors for our set of light and heavy quark masses, from the analytic method 
and using the conserved vector current. 

We are interested in the form factors at different q 2 , for quark masses corre- 
sponding to the physical D and K (ir) mesons. Thus we have to extrapolate the 
form factors reported in tables ^ and |], both in mass and momentum. We have 
proceeded as follows: 

i) At fixed heavy quark mass and light meson momentum, Pk,k*, the generic 
form factor F (F = f + , A\ ,...) has been extrapolated linearly in the light quark 
mass to values corresponding to the strange (D — > K, K*) or massless (D — > rr, p) 
quarks. In ref.Q, it was shown that SU(3) symmetry breaking effects are very 
small, i.e. that a linear dependence of the type: 



1.2 

q 



describes well the behaviour of F as a function of the light quark mass (m 
1/2 a x (1/K 12 — 1/K cr )). In eq.(|7|) K\ and K 2 correspond to the masses of the 
final and spectator quarks respectively. Thus the form factors extrapolated to 
the strange and light quark masses, F K and F n , read F K = a + (3(1/K s + 1/K cr ) 

4 We recall that all our calculations have been done with pjj = 0. 
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K U K 2 


V 


q 2 a 2 


f + (q 2 ) 


/v) 


0.1275 

U. ±^±00 


0,0,0 
1 n n 

1,1,0 


0.296(4) 
-0.010(2) 


U.Ol 1 O 1 

0.76(17) 


0.91(8) 
f) 74f7'l 

0.77(17) 


0.1275 

fl 14QD 


0. 0.0 

1 n n 

1, u, u 

1,1,0 


0.325(5) 
-0.003(3) 


fl 89H 1 
0.82(21) 


0.90(9) 
0.83(21) 


0.1275 
n 14Q^ 


0. 0.0 

1 n n 

1, u, u 

1,1,0 


0.361(7) 
n i ^1 (a\ 

0.005(3) 


fl 89H 7\ 

0.97(26) 


0.91(11) 
D 71 fl 1\ 

0.97(26) 


0.1325 

U. ±^±00 


0,0,0 
1 n n 

1,1,0 


0.186(3) 

U.UO I 1 i 1 

-0.078(1) 


fl 7Q(7\ 
0.74(17) 


0.93(8) 
n 7(\(7\ 

0.80(16) 


0.1325 
D 14QD 


0. 0.0 
1 n n 

1, u, u 
1,1,0 


0.208(4) 
-0.074(2) 


D 7R(1 r\\ 

U . I Ol 1U 1 

0.79(21) 


0.92(10) 
0.87(21) 


0.1325 
n 14Q^ 


0,0,0 
1 n n 

-L , U, U 

1,1,0 


0.236(6) 
-0.069(2) 


n 7%('\ ^\ 
0.94(26) 


0.93(11) 
1.03(28) 


0.1375 
n 148^ 


0,0,0 
1 n n 

1,1,0 


0.097(2) 
D D97M 

-0.122(1) 


D 7^7"! 
u. t oy i j 

0.70(17) 


0.95(9) 
0.83(16) 


0.1375 

D 14QD 


0,0,0 
1 n n 

1,1,0 


0.112(3) 
n D91 M \ 

-0.120(1) 


n 7'Uq'i 

U. t O\o J 

0.75(21) 


0.95(10) 
0.92(22) 


0.1375 


0,0,0 
1 n n 

1,1,0 


0.132(4) 
-0.118(1) 


D 79M "V\ 
0.89(26) 


0.95(11) 
D 74M4^ 

1.10(31) 


0.1425 
fl 1485 


0,0,0 

1 n n 

1,1,0 


0.033(1) 
-fl C\f,1(-) 

U.UUtJ I J 

-0.137(-) 


0.65(17) 


0.97(9) 
0.88(17) 


0.1425 
0.1490 


0,0,0 
1,0,0 
1,1,0 


0.042(1) 
-0.061(-) 
-0.137(-) 


0.66(8) 
0.68(21) 


0.96(10) 
0.77(10) 
0.98(25) 


0.1425 
0.1495 


0,0,0 
1,0,0 
1,1,0 


0.053(2) 
-0.059(1) 
-0.137(-) 


0.62(11) 
0.80(27) 


0.96(11) 
0.75(15) 
1.16(40) 



Table 3: Pseudoscalar — > pseudoscalar form factors for different momenta and 
quark masses using the "analytic" method. The values reported in the table have 
been obtained using the "conserved" current. Errors denoted by (-) are smaller 
than the figures reported in the table. We also give the squared momentum transfer 
in lattice units. 
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K U K 2 


P 


q 2 a 2 


V(q 2 ) 


A(g 2 ) 


A 2 (q 2 ) 


A(q 2 ) 


0.1275, 
0.1485 


0,0,0 
1,0,0 
1,1,0 


0.237(6) 
0.088(4) 
-0.035(3) 


0.80(7) 
0.61(15) 


0.61(6) 
0.52(4) 
0.47(15) 


0.43(12) 
0.36(34) 


0.09(3) 
-0.04(3) 


0.1275, 
0.1490 


0,0,0 
1,0,0 
1,1,0 


0.254(8) 
0.096(5) 
-0.030(3) 


0.77(10) 
0.59(18) 


0.59(6) 
0.51(5) 
0.44(20) 


0.45(16) 
0.30(42) 


0.12(5) 
-0.03(3) 


0.1275, 
0.1495 


0,0,0 
1,0,0 
1,1,0 


0.272(9) 
0.105(5) 
-0.025(4) 


0.76(16) 
0.56(24) 


0.56(7) 
0.52(8) 
0.37(28) 


0.51(26) 
0.19(57) 


0.17(8) 
-0.03(4) 


0.1325, 
0.1485 


0,0,0 
1,0,0 
1,1,0 


0.139(5) 
0.011(3) 
-0.095(2) 


0.82(7) 
0.62(16) 


0.66(6) 
0.56(5) 
0.50(16) 


0.43(12) 
0.34(31) 


0.013(7) 
-0.12(9) 


0.1325, 
0.1490 


0,0,0 
1,0,0 
1,1,0 


0.152(5) 
0.016(3) 
-0.092(2) 


0.79(10) 
0.58(19) 


0.64(7) 
0.55(6) 
0.47(21) 


0.44(16) 
0.29(39) 


0.02(1) 
-0.12(12) 


0.1325, 
0.1495 


0,0,0 
1,0,0 
1,1,0 


0.165(7) 
0.022(4) 
-0.089(2) 


0.78(17) 
0.55(24) 


0.60(8) 
0.56(9) 
0.39(31) 


0.48(25) 
0.16(53) 


0.04(2) 
-0.11(16) 


0.1375, 
0.1485 


0,0,0 
1,0,0 
1,1,0 


0.064(3) 
-0.042(1) 
-0.130(6) 


0.84(8) 
0.62(17) 


0.71(7) 
0.59(5) 
0.53(16) 


0.43(10) 
0.32(27) 


-0.06(2) 
-0.19(14) 


0.1375, 
0.1490 


0,0,0 
1,0,0 
1,1,0 


0.072(4) 
-0.039(2) 
-0.129(1) 


0.81(11) 
0.57(20) 


0.69(7) 
0.58(7) 
0.50(22) 


0.43(15) 
0.26(35) 


-0.07(3) 
-0.19(19) 


0.1375, 
0.1495 


0,0,0 
1,0,0 
1,1,0 


0.081(5) 
-0.036(2) 
-0.128(1) 


0.79(17) 
0.51(24) 


0.64(8) 
0.59(10) 
0.41(32) 


0.43(23) 
0.10(48) 


-0.08(4) 
-0.18(27) 


0.1425, 
0.1485 


0,0,0 
1,0,0 
1,1,0 


0.015(1) 
-0.068(2) 
-0.135(-) 


0.87(9) 
0.61(18) 


0.75(8) 
0.62(6) 
0.56(17) 


0.41(8) 
0.27(22) 


-0.13(5) 
-0.26(18) 


0.1425, 
0.1490 


0,0,0 
1,0,0 
1,1,0 


0.019(2) 
-0.067(-) 
-0.136(-) 


0.81(12) 
0.53(21) 


0.72(8) 
0.61(8) 
0.52(22) 


0.39(12) 
0.20(28) 


-0.15(6) 
-0.26(24) 


0.1425, 
0.1495 


0,0,0 
1,0,0 
1,1,0 


0.023(2) 
-0.066(1) 
-0.136(-) 


0.78(18) 
0.43(25) 


0.67(9) 
0.60(11) 
0.40(33) 


0.36(19) 
0.02(40) 


-0.17(10) 
-0.21(35) 



Table 4: Pseudoscalar — > vector form factors for different momenta and quark 
masses using the "analytic" method. The values reported in the table have been 
obtained using the "conserved" current. We also give the squared momentum 
transfer in lattice units. 



9 





K 2 


aM 5 


aM v 


aM A 


aM s 


0.1275 


0.1485 


0.813(4) 


0.822(4) 


0.94(2) 


0.91(2) 


0.1275 


0.1490 


0.804(4) 


0.812(4) 


0.93(3) 


0.90(3) 


0.1275 


0.1495 


0.794(5) 


0.803(5) 


0.92(3) 


0.90(4) 


0.1275 


K s 


0.795(6) 


0.803(6) 


0.92(3) 


0.90(4) 


0.1275 


K 


0.773(6) 


0.783(5) 


0.91(5) 


0.88(6) 


0.1325 


0.1485 


0.700(4) 


0.712(4) 


0.83(2) 


0.81(3) 


0.1325 


0.1490 


0.689(4) 


0.702(4) 


0.83(3) 


0.80(3) 


0.1325 


0.1495 


0.679(5) 


0.693(5) 


0.82(4) 


0.79(5) 


0.1325 


K s 


0.680(6) 


0.693(6) 


0.82(4) 


0.79(5) 


0.1325 


K 


0.657(6) 


0.672(6) 


0.81(5) 


0.77(7) 


0.1375 


0.1485 


0.580(4) 


0.598(5) 


0.73(3) 


0.69(3) 


0.1375 


0.1490 


0.569(4) 


0.588(5) 


0.72(3) 


0.68(4) 


0.1375 


0.1495 


0.558(5) 


0.578(5) 


0.72(4) 


0.67(6) 


0.1375 


K s 


0.558(6) 


0.578(7) 


0.72(4) 


0.67(5) 


0.1375 


K cr 


0.533(6) 


0.556(6) 


0.71(6) 


0.65(8) 


0.1425 


0.1485 


0.450(4) 


0.479(5) 


0.63(3) 


0.56(4) 


0.1425 


0.1490 


0.438(4) 


0.468(6) 


0.62(4) 


0.55(5) 


0.1425 


0.1495 


0.425(5) 


0.457(6) 


0.62(5) 


0.54(7) 


0.1425 


K s 


0.425(6) 


0.457(8) 


0.62(5) 


0.54(7) 


0.1425 


K cr 


0.396(5) 


0.433(7) 


0.60(7) 


0.51(10) 



Table 5: Pseudoscalar (M 5 ), vector (My), axial (Ma) and scalar (M s ) meson 
masses for mesons composed of heavy-light quarks (dimensionless units). 
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and F n = a + 2/3/ K cr , where K s = 0.1495(1) is the value of the Wilson parameter 
for the strange quark and K cr = 0.1506(1) the critical one. The value of K c h arm 
corresponding to the charmed quark is 0.1379(11)[|. 

ii) Next we have extrapolated F K,7T in the mass of the heavy quark according 
to the expression: 

F**-A + JL, (8) 

where Mp is the mass of the heavy meson. At each value of Pk,k* , we have also 
extrapolated in 1/Mp to the D and B mesons using the dependence expected in 
the heavy quark effective theory (HQET) |T7| , see eq.flTTD below. The difference of 
the results obtained with the two methods are discussed later on. For D decays, 
we have verified that extrapolating in a different order, e.g. first in q 2 , then in 
the light quark mass, and finally in the heavy quark mass, leads to very similar 
results, within the statistical errors. 

iii) To obtain the form factors at q 2 = 0, we have extrapolated the form 
factors at pk,k* = (1, 0, 0) (which in many cases corresponds to the smallest q 2 , 
cf. tables || and ^) by assuming meson dominance: 



F(0) 
q 2 /M 2 



where F can be / + , V, A\ etc. and M t is the mass of the lightest meson exchanged 
in the t-channel. Thus the vector, scalar and axial scalar meson masses have been 
used for the extrapolation of / + and V, f° and v4 12 respectively. M t is computed 
on the lattice, over the same configurations, at the same heavy and light quark 
masses used for the three-point functions. In table [5] we report, in lattice units, 
the heavy-light pseudoscalar, vector, scalar and axial masses used in the present 
analysis, at different values of the Wilson parameters. In table ^| we also report 
the physical masses as extrapolated from the values given in table [5| To the 
results in table | we add M Ds - M Dd = (93 ±6) MeV and M Bs -M Bd = (75 ± 12) 
MeV Q. In doing the extrapolation from small g 2 's to q 2 = 0, the precise value of 
M t is relatively unimportant. For example we have verified that by using in all 
cases the vector meson mass, the results change by about 5%. 

In order to show the q 2 dependence of the different form factors and to com- 
pare it with the meson dominance hypothesis, we display in figs.|I] the various form 

5 The errors given here are slighty different from those on ref. [g[ where the time interval 14-24 
was used for the fits of the two-point functions. In the present study we find a -1 = 3.6(2) GeV 
instead of 3.7(2) as in ref.^|. For the heavy-light mesons the interval 14 — 24 is probably more 
appropriate, but we have checked that the difference for the form factors is immaterial given 
the statistical errors of the results. 

6 In ref. Q were we used a time interval more apropriate to spectroscopy, we found Mr> s — 
M Dd = (86 ± 5) MeV and M Bs - M Bd = (62 ± 7) MeV. Using a improved action(l| in the 
static limit it has been found Mb„ — Ms d = (76 ± 10) MeV. These values can be compared to 



the experimental results M Ds - M Dd = (99.9 ± 0.7) MeV|19 and M Bs - M Bd ~ (96.3 ± 4. 



MeV obtained from M Bs = 5374.9 ± 4.4 MeV|2(| and M Bd = 5278.6 ± 2.0 MeVp 
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M D * 


M D ** 1++ 


Mb** 0++ 


AW 


M Da ** 1++ 


M Ds ** 0++ 


1.95 ±0.01 


2.50 ±0.22 


2.30 ±0.28 


2.04 ±0.02 


2.55 ±0.16 


2.38 ±0.19 


Mb* 


Mb** 1 ++ 


M B ** ++ 


M Bs * 


M Ba ** 1 ++ 


M Bs » 0++ 


5.27 ±0.01 


5.66 ±0.15 


5.61 ±0.17 


5.31 ±0.02 


5.72 ±0.10 


5.68 ±0.11 



Table 6: Masses in GeV predicted from the lattice for the vector, axial and scalar 
excitations of the D and B mesons (remember that the pseudoscalar masses Mb 
and Mb are used as an input). These masses have been used in several occasions, 
as mentioned in the text, to extrapolate the form factors at q 2 = 0. 



factors. They are shown as a function of the dimensionless variable q 2 /M 2 , for 
K H = 0.1375 and K L = 0.1495, which corresponds to the meson masses closest 
to the physical ones for D — > K, K* decays. For any given form factor, M t is the 
lattice meson mass appropriate for that particular channel. In the figures we have 
also included the points corresponding to a final meson with momenta (1,1,0) 
and (2, 0, 0). These points are reported for completeness but, given the extremely 
large systematic effects, they have never been used in the analysis. These figures 
show that the behaviour is compatible with meson dominance, at least in the 
range of masses and q 2 explored in our simulation. It is unknown whether the 
meson dominance will remain valid when we extrapolate to B decays, since the 
range of q 2 extends much further away from the pole in that case. Notice that 



in ref. [21 , using QCD sum rules, it was found that the axial form factors do not 



follow the behaviour expected on the basis of the meson dominance. 



2.3 Main Systematic Effects 

Before we proceed to the physics of B and D meson decays let us summarize and 
discuss the main sources of systematic effects, besides "quenching", present in 
our calculation: 

a) renormalization constants: the perturbative calculation is not unique in the 
sense that one has freedom in choosing the appropriate expansion parameter. It 
has been suggested that the bare coupling constant g% — 6//3 is not very suitable 
and that it should be replaced by an effective coupling g 2 ^ [0]-p4|. In our case, 
using the recipe of ref. |23], g 2 ff/g% turns out to be ~ 1.59, which leads, using 
the formulae of refs.[l4||-||16||, to Zy = 0.75 and Za = 0.84. The reliability of these 



estimates can be tested by the non-perturbative ratios (see also ref.[ 



<u\V£\P> 
<a\VM(3>> 



Z v = Z N7TT^ > (10) 



where V L ' are the local and conserved currents respectively. In our simulation, 
we found values of Zy which vary between 0.66 and 0.95 depending on the matrix 
element, see table [7]. This range of values reflects uncertainties which have been 
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Figure 1: We show an example of the q 2 behaviour of the form factors. We 
have chosen for the light and heavy Wilson parameters K=0.1J h 95 and 0.1375 
respectively. In each plot the two points furthest to the left (corresponding to 
the momentum assignment (1,1,1) and (2,0,0)) are only shown for the sake of 
illustration. We have never used them in the fits due to their large systematic 
and statistical errors. The curves correspond to the nearest pole dominance ap- 
proximation with the mass M t of the relevant meson taken from table [| and the 
numerator taken to fit the point closest to q 2 = i.e. the momentum assignment 
(1,0,0). 



interpreted as 0(a) effects |27| an d shown to be much smaller with the use of 
an "improved" action a la Symanzik||26||. In all our estimates we have used the 
perturbative values, Zy = 0.84 and Za = 0.88. With the present systematic 
uncertainties this choice is appropriate. The reader can easily replace our numbers 
by those corresponding to his preferred value for g 2 ff- 

b) 0(a) effects: from the above discussion we expect a systematic error of 
the order of 10-20% coming from 0(a) effects in the current matrix elements. 
The Fermilab group has suggested that this systematic error may be reduced by 
multiplying the propagators of the heavy quark Q by exp(m<ga) and redefining 
its mass[^8|. This prescription, which does not change the results for m^a <C 1, 
is motivated by the tree-level behaviour of lattice propagators. It is expected to 
work for small momenta, p^a <C 1, when the quark mass is large, mqa ~ 1. The 
difference between the results obtained using the standard procedure (i.e. multi- 
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p 


< ViVt > 


< KJ D > 


< KJiD > 


< K^JiD > 


(0,0,0) 


0.659(3) 


0.87(-) 






(1,0,0) 


0.652(3) 


0.94(1) 


0.75(1) 


0.66(2) 


(1,1,0) 


0.645(5) 


0.95(3) 


0.79(3) 


0.69(5) 



Table 7: Values of Zy calculated by taking the ratio of three-point correlation func- 
tions with inserted the "conserved" weak vector current divided by the correspond- 
ing correlation function with the local current inserted, cf. eq.(10). The Wilson 
parameters used for the three-point functions are Kh = 0.1375 and Kl = 0.1495. 
The results are shown for several momenta of the K, K* , and for different Lorentz 
components of the weak current, J i} and K* source, K* . We also give Zy as de- 
rived from the light-light two-point function in the column labelled as < V\V\ >, 
computed with Kl = 0.1495. 



plying the lattice quark propagators only by 2K), and those obtained following 
the modified prescription at least represent some measure of the uncertainty due 
to discretisation errors. 

We observe that the modification proposed by the Fermilab group leads to 
a universal change of all the matrix elements of a given current. Consequently, 
it cannot repair the fact that one finds very different values of Zy for different 
matrix elements as shown in table |7|. Furthermore we have checked that Zy from 
< ViVi > is almost independent of the quark mass. Indeed Zy from < ViVi > 
for heavy-light mesons ranges from 0.643(2) to 0.653(3), which is almost identical 
to the light-light case reported in table |7|. We conclude that in this case the 
dependence of Zy on the quark mass does not follow the behaviour predicted in 
ref. |[28|| []. For this reason we have not followed the suggestion of ref.[^8j. We 



have instead decided to report here the form factors obtained from the local 
and "conserved" vector current and the local axial current. At the end, we will 
average the results obtained from the local and conserved vector current and add 
(in quadrature) as a systematic error (~ 10 — 20% for V(0)) their difference. In 
the case of the axial current we only have results from the local axial current, for 
which we also expect a 10-20% systematic uncertainty coming from O(a) effects. 

c) extrapolation to the physical quark masses: the extrapolation in the light 
quark mass is quite natural and unlikely to be a source of an important uncer- 
tainty within our statistical accuracy. More delicate can be the extrapolation 
in the heavy quark mass. There are arguments, based on HQET which allow 
the expansion of the form factors at fixed Pk,k*, with \pk,k* \ *C Mp, in inverse 
powers of the heavy meson mass Mp[pjJ. All the relevant formulae are given in 



sec.EO, where we discuss the extrapolation to the B meson. Here we simply state 



7 This indicates that the vacuum to vector matrix element has a dependence on the quark 



mass different from that observed for the forward matrix element of the local current) 28 
Similarly one may worry about the use of the recipe of ref.pg]] for the pseudoscalar meson 
decay constants. 
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V{0) 


A 1 (0) 


A 2 (0) 


"a" -analytic 


0.73(72) ±0.16(16) 


0.85(1.10) ±0.24(30) 


0.63 ±0.14 


0.46 ±0.33 


"b" -analytic 


0.74(72) ±0.16(16) 


0.85(1.10) ±0.24(30) 


0.63 ±0.14 


0.46 ±0.33 


"a" -ratio 


0.60(59) ±0.12(12) 


0.84(1.09) ±0.24(30) 


0.63 ±0.14 


0.44 ±0.27 


"b"-ratio 


0.61(60) ±0.12(12) 


0.85(1.10) ±0.23(29) 


0.63 ±0.14 


0.44 ±0.27 



Table 8: Semi-leptonic form factors at zero momentum transfer for D — > K 
and K* using different extrapolations to the D mesons: eq.(8), labelled as "a", 
and eq.(ll), labelled as "b". The results are reported using the method called 
"analytic" or "ratio" to extract the form factors. The number obtained from the 
local vector current are reported in parenthesis. 



that HQET suggests that / + scales, at fixed Pk,k*, as: 

/+ = M^ 7+ x(l + ^) (11) 

— 1/2 

and similarly for /~, V and A2, while A\ scales like M P . For D mesons, the 
results, obtained by extrapolating the values given in tables ^| and [| by using 
eq.(|) (fit "a") or eq.(|ll]) ( fit "b") are reported in table § Q. Although we 
have found that the x 2 is slightly better if we fit the form factors according to 



eq. flTll) , we do not have sufficiently good data to distinguish between the two 
behaviours. For completeness we also give the results obtained by using the local 
vector current (in parenthesis). From this table one sees that for D mesons the 
difference between the extrapolation "a" and "b" is rather small, < 3%. There 
is more difference in the case of the vector current between the results obtained 
using the conserved or the local current, in particular in the case of D — > K*. 
As mentioned before we include this difference in the final error. For the axial 
current we are not able to estimate the error due to the determination of Za, 
since we do not have the equivalent of the conserved vector current. Besides the 
statistical errors reported in the tables, we then expect for the axial form factors 
a further error of order 10 — 20%. 



3 Physics Results 
3.1 D Meson Decays 

Given the discussion of the previous section, our best estimates for the form 
factors and partial widths are those reported in the abstract and in tables ^| 
and |l(J In this table we report our results together with other calculations 
and experimental determinations of the form factors. With respect to other 

8 Below we will give the same kind of table also for the B meson. 
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Ref. 


/£(o) 


V(0) 


Ai(0) 


A 2 (0) 


this work 


0.65 ±0.18 


0.95 ±0.34 


0.63 ±0.14 


0.45 ±0.33 


0-1 


1 


0.63 ±0.08 


0.86 ±0.10 


0.53 ±0.03 


0.19 ±0.21 


0-1 


1 


0.90 ±0.08 ±0.21 


1.43 ± 0.45 ± 0.49 


0.83 ±0.14 ±0.28 


0.59 ±0.14 ±0.24 


IP 


1 


0.76 


1.23 


0.88 


1.15 




I 


0.76-0.82 


1.1 


0.8 


0.8 


E 


I 


6 +tU5 


1.1 ±0.25 


0.5 ±0.15 


0.6 ±0.15 


Exp. |3| 


0.70 ±0.08 


0.9 ±0.3 ±0.1 


0.46 ±0.05 ±0.05 


0.0 ±0.2 ±0.1 


Ref. 


j\ i f (7 1 


V(0)Mi(0) 


A 2 (0)/A 1 (0) 


«/ V ^ mace / 


this work 


0.62 ±0.09 


1.50 ±0.28 


0.7 ±0.4 


0.93 ±0.13 


1 


0.77 ±0.20 


1.6 ±0.2 


0.4 ±0.4 




1 




1.27 ±0.16 ±0.31 


1.99 ±0.22 ±0.33 


0.7±0.16±0.17 






I 




1.1 


1.3 


1.15 


J31 


I 




1.4 


1.0 




m 


I 




2.2 ±0.2 


1.2 ±0.2 




Exp. g 


0.54 ±0.06 ±0.06 








Exp. g 




2.00 ±0.33 ±0.16 


0.82 ±0.23 ±0.11 





Table 9: Semi-leptonic form factors for D — > K and K* . For f^ and V we have 
averaged the results obtained from the conserved and local currents and considered 
the difference as an error to add in quadrature to the statistical one. 



calculations at a lower value of (3, for example refs. 0-01 the errors are larger. This 
is expected because the statistics is lower and f3 is higher. It is however reassuring 
that, at a value of a -1 almost twice that of refs.||-[|J, we find compatible results. 
Moreover experiments and lattice calculations are in good agreement for 
V and A\. On the other hand for A 2 the situation is unclear. E691 and ref.[|J 
suggest a smaller value of A 2 than E653, ref. and the present work, even though 
the errors on this quantity are so large that all the results are compatible. Thus 
for A 2 it is very important to reduce both the experimental and theoretical errors. 
The agreement with the other form factors is already remarkable, given the fact 
that the lattice calculations have no free parameter. 

3.2 Extrapolation to B meson decays 

With the present lattice spacing we are unable to study directly the B mesons 
in numerical simulations. However indirect information on B physics may be 
available using the following strategy. One first studies the scaling behaviour 
of a given physical quantity in the region of the charm quark mass. Then one 
computes the same quantity in the effective ("static") theory, i.e. in the limit in 
which the heavy quark mass is infinite. By interpolating between the two results 
with the help of the expected scaling laws, one can predict that physical quantity 
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Ref. 


T(D ^ K)/10 10 s- 1 


T(D -> iT^/lO 10 ^ 1 


rp^^/io 10 ^- 1 


this work 


6.8 ±3.4 


6.0 ±2.2 


0.56 ±0.36 




1 




5.8 ±0.15 


5.0 ±0.9 


0.5 ±0.2 


exp. 


Pi 


3,|32| 


7.0 ±0.8 


4.0 ±0.7 




Ref. 


T(D -> p)/ltf a 8- x 


r(D -> K*)/T(D -> X) 


ryr r 


this work 


0.50 ±0.23 


0.92 ±0.55 


1.27 ±0.29 


I 




0.4 ±0.09 


0.86 ±0.22 


1.51 ±0.27 







1.14 


0.89 







1.45 


1.11 


m 




0.5 ±0.15 


0.86 ±0.06 


® 




0.63 ±0.09 




exp. 


35| 




0.57 ±0.08 


1.15 ±0.17 



Table 10: Semi-leptonic partial widths for D —> K , K* , it and p, using V cs = 
0.975 and V c d = 0.222. We also report the ratio of the longitudinal to transverse 
polarisation partial widths for D — ► K * . 



for the B meson. The value in the static limit reduces the uncertainty due to the 
extrapolation from the charm region. This strategy has proved to be effective for 
the pseudoscalar decay constant, see for example refs.|| |2"9"| . 

At our large value of j3, the small statistics is not sufficient to compute semi- 
leptonic form factors in the static limit. Nevertheless we can study the scaling 
behaviour of the form factors and try an extrapolation to the b quark. With 
the actual errors, our predictions remain necessarily at a semi- quantitative level. 
The study is however interesting in itself. We have found that the expected 
dependence of the form factors in l/M P , see eq . , is compatible with our 
results. For some of the form factors however the corrections to the static results 
in the charm region may be sizeable. In the following we explain the extrapolation 
of the form factors, computed at several values of the heavy quark mass, to the 
B meson. 

On the basis of HQET, up to 0(1/ Mp), up to logarithmic corrections, one 



expects the following behaviour for the relevant form factors ||17|| : 

W = 7+x| [+ ^> W = " vx{+ ^ > 

/-, , ^2 \ . ,,1/2 (\ , $1 



The expansions given in eqs. (|I~2D become valid in the limit of large rriQ, at fixed 
momentum p of the light meson (in the frame where the heavy meson is at rest) 
and when \p\ <C mg ~ Mp. The above conditions are always satisfied for q ma xi 
when the initial and final mesons are both at rest. Provided we are not too close 
to the meson pole in the t-channel, the dependence on the light mass is smooth 
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Figure 2: Form factors extrapolated to the chiral limit for the light quark, as a 
function of the inverse pseudoscalar mass (1/Mp ) for the momentum assignment 
(1,0,0). The crosses are the lattice points, the diamonds are the extrapolation to 
the D and B meson. Notice that the D is very close to a lattice point. The points 
corresponding to the lightest heavy quark mass (furthest to the right) has not been 
used in the fits. 



and it is well described by eq.(0). We can then extrapolate in the heavy quark 
mass according to eqs.flT2|). The points corresponding to p = (1, 0, 0) also satisfy 
the above conditions. However, in the range of masses where we compute the 
form factors, niq ~ m c ^ arm , these points correspond to q 2 ~ 0. This is due to the 
fact that the heavy meson masses are not so heavy in comparison with \p\ and 
that the light mesons are not really light enough. Thus, while the typical yfq 1 on 
our lattice is at most 1.4 GeV, the extrapolation in the heavy quark mass will 



bring us to yfq 1 = 4.2 GeV (q 2 = M% - 2 * M B ^/M 2 + (2vr/24a) 2 + M 2 ), and 
similarly for B — > p. The uncertainty involved in such an extreme extrapolation 
must be borne in mind. The validity of the extrapolation is partially justified only 
because the range in 1/Mp is relatively small and most of the form factors have 
a smooth behaviour. The case (1,1,0) is more complicated because, with the 
heavy quark masses at hand, we have \p\ ~ Mp where the expansion of ref.[17] is 
questionable. Anyhow the errors are so large that we have not used these points. 
In fig.^| we give the form factors, at the value of q 2 corresponding to the 
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p 


7+ GeV" 1 / 2 


7y GeV" 1/2 


7i GeV +1 / 2 


7 2 GeV~ 1/2 


(0,0,0) 






0.96 ±0.16 




(1,0,0) 


0.39 ±0.25 


0.29 ±0.12 


1.05 ±0.25 


0.44 ±0.25 


p 


5 + GeV 


5 V GeV 


5 X GeV 


5 2 GeV 


(0,0,0) 






-0.33 ±0.09 




(1,0,0) 


0.0 ± 1.1 


1.9 ± 1.3 


-0.46 ±0.22 


-0.6 ±0.8 



Table 11: The coefficients of the 1/rriQ expansion of the form factors defined in 
eqs.(12). 



momentum assignment (1, 0, 0), extrapolated to the chiral limit in the light quark 
mass, as a function of 1/Mp (crosses). In the figure we also give the value 
extrapolated to the D and B mesons (diamonds). The figure shows that the 
points are smoothly extrapolated, even in the case of V since 7y is small. This is 
encouraging for further studies with higher statistics. In table |ll| we report the 
values of 7+,v;i,2 and £+, 1/1,2 for p = (0,0,0) and (1,0,0), in physical units. We 
tried also (1, 1,0), but the errors turned out be of order 100% and we have not 
reported them in the table. We notice that the first correction in 1/Mp is small 
in the case A\. This is also the case for / + , at (1,0,0), even though, because of 
the large error, we cannot exclude a slope ~ 0.7 GeV. The same can be said for 
A2, whose value anyhow is badly determined even before the extrapolation. In 
the case of V we find instead a rather large correction. We believe that this is 
an interesting and rather unexpected result, which deserves an effort to reduce 
the size of the statistical error. Taking our present errors into account we cannot 
exclude that the large value of 5y is correlated to the small value of 7y and that 
it will disappear with more accurate results. 

From the numbers given in the table we can predict the form factors of the B 
mesons. We give all the form factors at q 2 = 0. They have been obtained from 
the form factors at p — (1,0,0), by using the meson dominance, with the mass 
of the mesons exchanged in the t-channel reported in table || These masses have 
been obtained by fitting the mass difference AM = Mp* — Mp (vector case) and 
AM = Mp** —Mp (axial case) as AM = Am + Bm/Mp. The results are reported 
in table |6|. They are rather close to the experimental masses when known, and 
also to the masses computed in ref.flBTl, i.e. M B * = 5.32 GeV and M B ** x ++ = 5.71 



GeV, except for the scalar mass which in ref.[|30| is very large: Mp** t0 ++ = 5.99 
GeV. The differences between the extrapolations using different sets of masses 
turn out to be small and may be smaller than the error induced by the assumption 
of the meson dominance. The results are reported in table 0, label "b", together 



with the results of refs.[|3"0"|-[]3"2" for comparison. In the evaluation of the errors 
we have also taken into account of the difference between the results obtained 
using the local or the conserved current. To show the stability of the results with 
respect to a different extrapolation, we have also reported the values obtained 
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Ref. 


/ + (0) 


V(0) 


A x {0) 


A 2 (0) 


this work "a" 


0.28 ±0.14 


0.37 ±0.14 


0.24 ±0.06 


0.39 ±0.24 


this work "b" 


0.33 ±0.17 


0.40 ±0.16 


0.21 ±0.05 


0.47 ±0.28 




33 


0.33 


0.33 


0.28 


0.28 


[31] 


0.09 


0.27 


0.05 


0.02 






0.26 ±0.02 


0.6 ±0.2 


0.5 ±0.1 


0.4 ±0.2 


Ref. 




V(0)/A 1 (0) 


A 2 (0)/A 1 (0) 




this work "a" 




1.4 ±0.5 


1.5 ±0.8 




this work "b" 




1.7 ±0.6 


2.3 ± 1.1 




\m 




1.0 


1.0 





Table 12: Semi-leptonic form factors for B —> tx and p. For f + and V we have 
averaged the results obtained from the conserved and local currents and considered 
the difference as an error to add in quadrature to the statistical one. The label "a" 
refers to the naive extrapolation in 1/Mp, eq.(8), and label "b" to the extrapolation 
given in eqs.(12) respectively. To extrapolate to zero momentum transfer we have 
used the masses of table 6. 



with the naive scaling given in eq. (§) , labelled as "a" . 

In the D meson the difference between the HQET scaling laws and naive 
scaling was immaterial. The differences remain small in the B case and with 
our errors we cannot distinguish the two behaviours. Notice that, because of the 
different scaling laws, A 2 j A\ increases with the mass of the heavy quark and can 
be ~ 2 for B — ► p. 



For most of the form factors the predictions of ref. [31| are much lower than 



all the others, cf. table |12[ These result in a much larger estimate of \V u b\, for 
a given experimental branching ratio. Notice however that in ref. [3~I] the form 
factor is computed at q^ax and then a "tempered" exponential dependence on 
q 2 is assumed. This q 2 yields a dramatic suppression at small q 2 for B meson 
decays, where the range in q 2 is very large. 

From the numbers reported in table [T2|, we can give a very rough estimate of 



the B — > 7r and B — > p branching ratios. The branching ratios are obtained from 
the form factors in the table, by assuming meson dominance for their dependence 
on q 2 , see eq.(^j). To get an estimate of the errors, we have allowed the form factors 
to vary in all the possible ways by one a within the statistical errors and to vary 
in all possible ways among the values obtained with different extrapolations in 
1/Mp, fits "a" and "b". In this way we can partially account for the uncertainty 
coming from the extrapolation and the q 2 dependence of the form factors. We 
finally get: 

T(B -> -kIvi) = |Kb| 2 (12 ± 8) x lO 12 ^ 1 (13) 

and 

T(B -> pLv{) = |Kfe| 2 (13 ± 12) x lO 1 ^- 1 (14) 
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Thus the errors are still too large to give more than an upper bound. Since 
experiment, on its own side, only gives upper bounds for the considered branching 
ratios, we cannot with our present accuracy extract any information on V u b- We 
mainly want to stress the feasibility of this analysis, provided one disposes of 
larger statistics. 
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